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We numerically investigate the role of quantum fluctuations in superresolution of 
optical objects. First, we confirm that when quantum fluctuations are not taken 
into account, one can easily improve the resolution by one order of magnitude 
beyond the diffraction limit. Then we investigate the standard quantum limit of 
superresolution which is achieved for illumination of an object by a light wave in 
a coherent state. We demonstrate that this limit can be beyond the diffraction 
limit. Finally, we show that further improvement of superresolution beyond the 
standard quantum limit is possible using the object illumination by a multimode 
squeezed light. 



Last years have witnessed an increasing; interest to investigation of quan- 
tum effects in optical imaging ^21, One of the questions recently reconsidered 
in the light of this latest development is about the ultimate quantum limits 
of resolution in optical systems. A classical resolution criterion formulated 
at the end of the last century by Abbe and Rayleigh states that the optical 
resolution is limited by diffraction present in any optical system due to the 
wave nature of light. This diffraction limit was introduced by Rayleigh for 
a simple observation of a diffracted image by a human eye. However, nowa- 
days using the modern CCD cameras for detection of optical images with 
subsequent electronic treatment of the digitized signals one can often improve 
the resolution beyond the limit imposed by diffraction. Such superresolution 
techniques use some a priori information about the input object and are lim- 
ited not by diffraction but by different kinds of noise in the detection and 
the electronic reconstruction systems. It was recently demonstrated ^ that 
the ultimate limit of superresolution is determined by quantum fluctuations 
of light, and that the use of special kind of spatially multimode squeezed light 
should allow to increase the capability of superresolution schemes. 

In this paper we numerically simulate the role of quantum fluctuations for 
superresolution of two simple optical objects placed close to each other so that 
they cannot be resolved according to the Rayleigh criterion. We consider a 
simple one-dimensional scheme of diffraction-limited coherent optical imaging 
shown in Fig. The object of finite size X is situated in the object plane. 
The first lens Li performs the Fourier transform of the object into the Fourier 
plane where a pupil of finite size d is located. Diffraction on this pupil is a 
physical origin of the finite resolution distance in the scheme. The second lens 
L2 performs the inverse Fourier transform and creates a diffraction-limited 
image in the image plane. As mentioned above, to achieve superresolution 
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Figure 1. Optical scheme of one-dimensioned diffraction-limited optical imaging. 



one needs some a priori information about the object. In our case we know a 
priori that the object is confined within the area of size X and is identically 
zero outside. The spatial Fourier transform of such an object is an entire 
analytical function. Therefore, knowing the part of the Fourier spectrum 
within the area d of the pupil allows for an analytical continuation of the total 
spectrum and, therefore, for unlimited resolution. However, this analytical 
continuation is extremely sensitive to the noise in the diffracted image and, 
as will be illustrated below, is limited by the quantum fluctuations of light. 

To simplify notations we shall use the dimcnsionless coordinates s — 2x /X 
in the object and the image planes and the dimcnsionless coordinate ^ = 
2y/d in the pupil (Fourier) plane. Let the classical dimcnsionless complex 
amplitudes of the electromagnetic field in the object, Fourier, and the image 
planes be respectively a(s), and e(s). The amplitudes a(s) and /(^) are 

related by the Fourier transform performed by the lens ii. 



(1) 

where c = ndX/ (2Xf) is the space-bandwidth product of the imaging system. 

The second lens L2 performs the inverse Fourier transform and creates an 
image in the image plane. The object and the image complex amplitudes are 
related by an integral operator, 

, , sinfcfs — s')l / /N , / /„v 

eis)^J_^^±-Jiais')ds'. (2) 

The orthonormal eigenfunctions of this operator are given by 

ipkis) = -^Vfc(s), \s\ < 1, (3) 



where ipk (■5) are the prolate spheroidal functions and Xk are the corresponding 
eigenvalues^. To achieve superresolution one can decompose the object field 
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in the basis of the eigenfunctions c^fc(s) as 

oo 

a{s) ^^ak(pk{s), |s| < 1, (4) 
with the coefBcients calculated by 

Ofc^y a{s)ipk{s)ds. (5) 
Similar decomposition can be written in the image plane, 

oc 

e(s) = ^ Ckil^kis), -oo < s <(X), (6) 

fc=0 

and in the Fourier plane, 

oo 
fc=0 

The coefficients Ck are calculated as 



Cfc = / e(s)V'fc(s)ds, (8) 

J — OO 

while the coefficients fk are given by 

fk - /' /(Ov'fcCOrfe (9) 
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Using the properties of the prolate functions ^pk it can be shown that the 
coefficients Ck and fk are expressed through Uk as follows: 



Bk = V^kO-k, (10) 



/fc = iVAfcafc. (11) 

Therefore, detecting the image e(s) in the image plane using, for example, 
a sensitive CCD camera, one can calculate the coefficients Ck according to 
© and than reconstruct exactly the coefficients Uk of the object using H1U|I . 
Alternatively, one can set the CCD camera in the Fourier plane to detect /{£,), 
evaluate the coefficients fk according to (|5J and reconstruct ak using (|lf |l . The 
first method could be called superresolving microscopy, while the second one 
superresolving Fourier-microscopy since one detects the Fourier spectrum. It 
should be noted that, since in both cases we need the complex field amplitudes 
and not the intensities, one should use the homodyne detection scheme with 
a local oscillator. 

In our numerical simulations we have tried both detection schemes and 
have given preference to the Fourier-microscopy since it involves integration 
over the finite region of the pupil, while detecting the image e(s) requires 
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integration over an infinite area in the image plane. It turns out that due to 
oscillating behavior of the prolate functions ipkis) one needs to take unrealis- 
tically large area in the image plane to achieve significant superresolution. 

For numerical simulations we have taken a simple object of two Gaussian 
peaks, 



a{s) = exp 



2a2 



exp 



2a2 



< 1, 



(12) 



of width a separated by distance 2so- We choose 2so = 1 and ct = 0.1, so that 
two peaks are well separated in the input object. The Rayleigh resolution 
distance R — nX/{2c) in dimensionless coordinates is equal to tt/c, where c 
is the space-bandwidth product. In our simulations we work with c = 1. In 
this situation for 2so < we are beyond the Rayleigh limit. This is clearly 
seen in Fig. [3 where we have shown the input object and its image observed 
in the image plane. 



(a) 



(b) 




Figure 2. Double-peak object a{s) used in numerical simulations (a) and its image e{s) (b). 



Therefore, observing the image in Fig. |2J3, it is impossible according to 
the Rayleigh criterion to resolve two Gaussian peaks in input object. How- 
ever, applying the superresolution technique one can easily reconstruct the 
input object beyond the diffraction limit. Wc illustrate the result of such a 
reconstruction in Fig. |31 In this figure we show the exact Fourier spectrum 
of the input object, drown by a solid line, as a function of dimensionless co- 
ordinate f in the Fourier plane. Only a part of this spectrum shown by a 
bold line, within the area of the pupil, |^| < 1, is transmitted to the image 
plane. This is a reason of very large diffraction spread in the image plane 
shown in Fig. IJb. Three dotted bold lines in Fig. |3 correspond to the Fourier 
spectrum of the reconstructed object using 2, 4, and 6 prolate functions. We 
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c;an sec; that the rcicoiistructcd spectrum approaches the exact one; for ever 
higher spatial frequencies |^| as the number of prolate functions increases. 
With 6 prolate functions two spectra are very close to each other for spatial 
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Figure 3. Spatial Fourier spectrum of the object from Fig. 2a (solid line), its part trans- 
mitted trough the pupil (bold solid line), and the spectra reconstructed with 2, 4, and 6 
prolate functions (three dotted bold lines). 

frequencies |^| < 8. This corresponds to a superresolution factor of 8 over the 
Rayleigh limit. 

Up to now we did not take into account the fluctuations in the detection 
of the Fourier components by a CCD camera in the Fourier plane. However, 
such fluctuations are always present in the detection scheme due to technical 
imperfections and the quantum nature of light. The quantum fluctuations 
of light set the ultimate limit of superresolution in optical imaging. The 
quantum theory of the optical imaging scheme in Fig. 1 was developed in 
Ref. [3]. Applying the same methods for the Fourier- microscopy, we can write 
the photon annihilation operators a(s) and /(^) in the object and the Fourier 
plane as 

CxD OO 

K^) =^Y^(^k^>k{s) ^^^hkXk{.s), (13) 
fe=o fe=o 

oo oo 

hi) ^J2hM0 + J29kXk{0- (14) 

fe=0 fe=0 

Here Xfc are the orthonormal basis functions in the region |s| > 1 and 
1^1 > 1 introduced in Ref. [3], and bk and gk are the corresponding annihilation 
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operators. It can be shown that the operator- valued Fourier coefficients fk 
are given by 



fk = i'^iVXkak + Vl-Xkh)- 



(15) 



This relation is similar to the transformation performed by a beam-splitter 
with amplitude transmission coefficients i'^\/Xk and reflection coefficients 
i'^Vl ~ ^k, and preserves the commutation relation of the annihilation and 
creation operators in the Fourier plane. 

We can use Eq. H15|l for calculation of the coefficients d^^^ in the recon- 
structed object as 



fk 



'At 



1-A, 



Xk 



(16) 



where the superscript (r) indicates "reconstructed". As follows from Eq. 
the reconstruction of the input object is no longer exact because of the second 
term in Eq. 1)16(1 . This term contains the annihilation operators bk responsible 
for the vacuum fluctuations of the electromagnetic field in the area outside 
the object. It is important to notice that these vacuum fluctuations prevent 
from reconstruction of the higher and higher coefficients ak in the object 
because of the multiplicative factor ^(1 — Xk)/Xk- Indeed, the eigenvalues 
Xk become rapidly very small after the index k has attained some critical 
value. This leads to rapid "amplification" of the vacuum fluctuations in the 
reconstructed object that limits the number of the reconstructed coefficients 
Uk ■ Below we illustrate numerically the role of these quantum fluctuations in 
superresolution. 




Figure 4. Reconstruction of the spatial Fourier spectrum of the object with Ught in a 
coherent state with mean photon number (N) = 10^^. Five dotted Unes correspond to five 
random Gaussian realizations of the quantum fluctuations. 
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The relative value of quantum fluctuations depends on the signal-to-noise 
ratio in the input object which for the light in a coherent state is determined 
by the total mean number of photons passed through the object area during 
the observation time. For example, for a laser beam with A 1064 nm and 
optical power of 1 mW, and observation time of 1 ms we obtain the mean 
photon number of (N) = 5.3 • 10""^^. 

In Fig. 0] we have shown the results of reconstruction of the spatial spec- 
trum of the object from Fig. |2K when quantum fluctuations of a coherent 
state are taken into account. The solid line gives an exact spatial Fourier 
spectrum of the object and the solid bold line the part of the spectrum passed 
through the pupil as in Fig. 13 We use 6 prolate functions and the mean 
photon number in the input object is taken (N) = 10^^. The five dotted lines 
correspond to five random Gaussian realizations of the quantum fluctuations 
in the coherent state of Ofe and the vacuum fluctuations of bk- The dotted bold 
line corresponds to the reconstructed spectrum by 6 prolate functions without 
noise (as in Fig. [2J|. One can observe that the role of quantum fluctuations 
becomes more and more important as one goes to the higher and higher spa- 
tial frequencies where the random realizations of the Fourier spectra deviate 
more and more from the mean value given by the dotted bold line. 

In Fig.Elwe have increased the total mean value of photons to (N) = 10^^. 
This corresponds to an increased signal-to-noise ratio in the input object and 
should allow for better superresolution. This is illustrated in Fig. |31 where we 
can reconstruct higher spatial frequencies in the Fourier spectrum as compared 
to Fig. H 
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Figure 5. Same as Fig. 4 but with (N) = lO^'"^. 

The same result can be achieved by using multimode squeezed light in- 
stead of increasing the power of the source illuminating the object. This is 
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illustrated in Fig. where we have used (A^) = 10^^ as in the Fig. 01 but con- 
sidered the light in a multimode squeezed state with the squeezing parameter 
= 10 instead of the coherent state. As the result the fluctuations in the 
higher spatial frequencies are decreased that gives better superresolution. 
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Figure 6. Reconstruction of the spatial Fourier spectrum of the object with Ught in mul- 
timode squeezed state with mean photon number (A'^) = 10^^ and the squeezing parameter 
e'' = 10. 



In conclusion, we have numerically investigated the role of quantum fluc- 
tuations in reconstruction of spatial spectra of optical objects. We have 
demonstrated that when quantum fluctuations are not taken into account 
one can achieve superresolution of about factor 10 over the Rayleigh limit. 
We have confirmed that the limit of superresolution is set by quantum fluc- 
tuations and depends on the signal-to-noise ratio in the input object. For a 
given signal-to-noise ratio one can further improve superresolution by using 
multimode squeezed light. 

This work was supported by the Network QUANTImEI (IST-200-26019) 
of the European Union. 
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